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Abstract:  This  paper  provides  insights  into  the  problems  of  designing  a  robot 

with  joint  and  link  flexibility.  The  relationship  between  the  deflection  of  the  robot 
under  gravity  is  correlated  with  the  fundamental  frequency  of  vibration.  We  consider 
different  types  of  link  geometry  and  evaluate  the  flexibihty  potential  of  different 
materials.  Some  general  conclusions  and  guidelines  for  constructing  a  flexible  robot 
are  given.  .  ... 
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1  Introduction 


For  many  years  people  have  been  interested  in  the  control  of  flexible  structures. 
Problems  caused  by  flexibility  show  up  frequently,  ranging  from  automated  machin¬ 
ery  that  has  to  pause  momentarily  to  allow  vibration  to  damp  out  to  instability 
in  non-colocated  control  loops.  Space-based  structures  and  robots  are  particularly 
susceptible  to  severe  flexibility  problems  because  they  tend  to  be  long,  slender  and 
extremely  lightweight. 

As  the  accurate  control  of  a  flexible  robot  is  a  complicated  task,  researchers  have 
begun  by  building  simplified  flexible  robots.  The  simplest  form  is  a  single  mass-beam 
system,  explored  by  many  researchers  including  [1,2].  Two  degree  of  freedom,  planar 
models  have  been  built  both  as  systems  with  stiff  links/flexible  joints  [3],  as  a  single 
flexible  link/fast  end  effector  [4]  system,  and  as  two  flexible  links  with  stiff  joints  [5]. 

Planar,  low  degree  of  freedom  robots  are  good  for  experimenting  with  control 
techniques  at  low  frequencies  of  vibration.  However,  they  cannot  reproduce  the 
complications  of  a  true  space-based  robot:  vibrations  with  coupled  bending  and 
torsional  modes  and  joints  exciting  modes  orthogonal  to  their  plane  of  action.  The 
practical  problems  of  controlling  such  an  arm  has  lead  researchers  to  building  three 
dimensional,  flexible  robots,  as  by  [6]. 

Unfortunately,  a  flexible  earth-based  robot  suffers  from  an  inseparable  problem: 
gravity.  Unless  you  build  in  some  form  of  compensation  to  eliminate  the  effects  of 
gravity,  your  robot  will  sag.  The  standard  way  to  get  around  sagging  is  to  make  the 
links  of  the  robot  directionally  stiff  with  a  higher  stiffness  vertically  than  horizontally. 
But  directional  stiffness  can  eliminate  some  of  the  very  effects  that  you  are  studying. 
Consider  a  standard  two  link  robot  model  with  three  degrees  of  freedom:  two  revolute 
joints  at  the  base  and  a  revolute  joint  at  the  shoulder  which  is  the  same  plane  as  one 
of  the  base  joints.  Rotation  of  the  base  or  movement  of  the  two  co-planar  joints  will 
excite  vibration  in  the  radial  sense,  but  only  rotation  of  the  base  will  excite  vibration 
perpendicular  to  this.  The  directional  stiffness  has  artificially  raised  the  natural 
frequency  in  a  direction  that  may  cause  the  most  problems  in  a  true  space-based 
arm. 

The  motivation  for  this  research  was  the  intent  to  design  and  build  an  anthropo¬ 
morphic,  two  link,  three  degree  of  freedom  test  fixture  that  would  exhibit  both  link 
and  joint  flexibility.  The  design  needed  to  behave  like  a  teleoperated  space  robot, 
complete  with  low  frequencies  of  vibration  and  coupled  modes.  This  paper  deals  with 
the  underlying  issues  that  came  up  while  designing  links  for  the  robot  that  could  give 
the  desired  flexibility  and  not  break.  The  first  issue  is  how  the  endpoint  deflection 
of  the  robot  under  gravity  relates  to  the  lowest  natural  frequency  of  vibration  of  the 
robot.  The  second  issue  deals  with  the  best  shapes  and  materials  to  be  used  for  the 
flexible  links  of  the  robot. 
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Figure  1:  Single  cantilever  beam  with  a  mass 


2  Vibration  Considerations 

When  building  a  robot  that  will  have  primary  modes  of  vibration  in  all  diiec-tlons, 
a  natural  design  issue  is  what  will  be  the  lowest  natural  frequency  of  the  robot.  A 
low  fundamental  frequency  has  several  advantages  and  one  clear  disadvantage.  The 
low  fundamental  frequency  is  easy  to  observe  and  record.  It  allows  the  higher  modes 
of  vibration  to  occur  at  frequencies  that  may  also  be  visible.  But  the  disadvantage 
is  that  a  flexible  arm  will  sag  under  gravity. 

It  turns  out  that  the  deflection  of  the  arm  ur. .  r  gravity  is  a  very  good  way  to 
estimate  its  natural  frequency  of  vibration,  and  vice-versa.  To  demonstrate  this,  we 
will  begin  by  deriving  the  relationship  between  the  endpoint  deflection  of  a  single 
beam  and  its  natural  frequency.  Then  we  demonstrate  that  for  real  two  link  flexible 
systems,  the  formula  relating  endpoint  deflection  to  natural  frequency  forms  a  useful 
estimate  of  the  system’s  natural  frequency. 


2.1  Single  Beam  Under  Gravity  Loading 

There  is  a  useful  relation  between  the  natural  frequency  of  vibration  and  the 
deflection  of  a  single  beam  urder  gravitational  loading,  mentioned  by  [7].  If  we 
consider  a  single  cantilever  beam  with  a  mass  (as  shown  in  Figure  1)  and  zissume 
that  it  behaves  as  a  Bernoulli- Euler  beam  we  have 


where  is  the  deflection  of  the  end  of  the  beam  and  the  mass  of  the  beam  is 
considered  to  be  negligible.  The  spring  constant  of  the  beam  that  relates  the  endpoint 
deflection  to  the  force  acting  at  the  end  can  be  written  as 


To  a  good  approximation,  the  lowest  natural  frequency  of  the  beam  is  given  by 

1  [k 


27rV  m 
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Figure  2:  Deflection  of  a  Single  Beam  under  Gravity 

where  fg  is  in  Hertz.  By  combining  equations  (2)  and  (3)  we  express  the  frequency 
of  vibration  of  the  cantilevered  beam  as 


(4) 


Hence  the  natural  frequency  of  vibration  of  the  beam  can  be  approximated  as  a 
function  of  its  deflection  under  gravitational  loading.  We  can  also  write  the  equation 
this  way: 


or  using  =  9.8  m/sec^,  we  have 


8t  = 


25 


47r2// 


Si  =  —  centimeters 

J  Q 


f9  = 


Hz 


(5) 


(6) 


(7) 


where  6i  is  in  centimeters  and  fg  is  in  Hertz.  This  is  displayed  graphically  in  fig¬ 
ure  2.  Note  particularly  that  frequencies  under  a  few  hertz  result  in  extremely  large 
deflections  under  gravity. 

2.2  Two  Beams  Under  Gravity  Loading 

Equations  (6)  and  (7)  are  useful  formulas  to  keep  in  mind  when  you  are  considering 
the  behavior  of  a  single  link  flexible  robot.  Now  consider  a  two  link  robot  modeled  as 
a  two  beam  system  with  a  joint  mass  m\  between  the  links  and  an  additional  payload 
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mass  m2  suspended  at  the  ends,  as  shown  in  Figure  3.  To  make  this  a  realistic  model 
of  a  robot  arm,  we  include  the  length  of  the  first  mass  as  Lj.  We  simplify  the  analysis 
with  two  assumptions.  We  assume  that  the  masses  of  the  beams  are  negligible  as 
compared  to  the  masses  of  the  joint  and  the  payload.  We  assume  that  the  beams 
behave  as  Bernoulli-Euler  beams  in  bending.  The  flexibility  matrix  of  the  system 
(see  [8])  is 


■^1  ■ 

Oil 

«12 

■  Fx  ■ 

.  <^2  . 

<221 

0.22 

.  ^2  . 

where  <5i,  82  are  the  deflections  at  mi  and  m2  from  the  forces  F\  and  F2.  Cij  refers  to 
the  deflection  at  mass  i  due  to  a  unit  force  at  mass  j.  The  values  of  a,j  can  be  found 
as 


1  3T,  3L^\ 

+  517  + 


1  31,  3 

0.2  -  02.  -  ;^[l  + 


(Z-i  -h  Lj){L2  +  Lj)] 


1  1  3(Z/2  T  Tj)  3  .  ^  r  \2i 


(9) 

(10) 

(11) 


where  we  have  substituted  Kx  =  3E2I1/ Lx^  and  K2  =  3E2I2I L2^ ■  Using  the  standard 
assumption  of  harmonic  motion  and  replacing  the  forces  Fx  and  F2  by  inertia  forces 
F,  =  —miSi  =  we  then  find  the  vibrational  frequencies  by  calculating  the 

determinant  and  solving  for  uj  from 


det 


(oiimi  - 

021^1 


0x21^2 

{022^2  —  ;7^) 


=  0 


which  can  be  solved  explicitly  for  /„  =  u!/2w  as 


(12) 


v/2 


fn  —  +  ^22^2  +  ((Oll^l  +  ^22^2)^ 


—  4(ana22n2im2  —  ax2^rnxm2))^^^'j  ^ 


(13) 
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I  he  endpoint  deflection  of  this  system  due  to  gravitational  loading  is 


8ttp  =  a2iTnig  +  022^2^? 

If  we  substitute  this  into  (4),  we  get  the  approximation 

1 

a2i"ii  +  o.22m2 


(14) 


(15) 


This  value  can  be  compared  to  the  expected  vibrational  frequency  from  equation  (13). 

In  fact,  it  is  easy  prove  that  fg  <  /„.  So,  assuming  that  the  vibrational  analysis 
is  reasonably  good,  then  the  quick  fg  calculation  will  be  a  lower  bound  for  the  actual 
lowest  frequency  of  vibration  of  the  system.  As  long  as  /„  is  not  too  much  greater 
than  fg,  fg  forms  a  useful  estimate  of  the  natural  frequency.  In  the  next  section  of 
the  paper,  we  will  empirically  demonstrate  that  fg  is  a  good  estimate  for  real  two 
beam,  two  mass  systems. 

Equations  (13)  and  (15)  depend  fundamentally  on  the  assumption  that  the  system 
can  be  considered  to  be  a  long  beam  with  a  mass  on  the  end.  In  fact,  when  m2  —*  oc 
01  Ki  — >  CO,  fg  ~  /„.  So  as  the  system  more  closely  resembles  either  a  single 
mass/beam  system  (where  mi  is  completely  negligible)  or  a  system  with  just  an  end 
mass  (where  A\  is  so  stiff  that  you  can  treat  K2  as  built  into  a  wall),  the  closer  it 
matches  the  ideal  case.  If  you  keep  reasonably  “balanced”  values  for  your  parameters, 
in  the  sense  that  each  beam  participates  in  the  vibration  and  neither  mass  strongly 
exceeds  the  other,  /„  ~  /j  is  a  good  approximation. 


2.3  Comparison  of  Vibrational  Formulas 

We  would  like  to  compare  the  values  of  fg  and  /„  to  determine  the  usefulness 
of  the  “tip  deflection  under  gravity”  approximation.  The  vibrational  frequency  /„ 
is  a  function  /„  =  /„(Li,  L2, Tj, mi, m2,  fci,  ^2)  which  is  too  complex  to  graph  easily. 
Instead,  we  pick  two  sets  of  parameters  and  see  how  varying  them  affects  the  ratio 
of  estimated  vibrational  frequency  to  actual  vibrational  frequency. 

Two  typical  cases  of  beam  configurations  are  displayed  in  Table  1,  one  aluminum 
and  one  steel.  These  values  were  chosen  as  representative  of  the  types  and  sizes 
of  systems  that  we  have  considered  in  the  course  of  our  research.  Each  of  these 
robot  configurations  has  a  natural  frequency  of  vibration  of  approximately  3.5  Hertz. 
Additionally,  the  K  values  were  choosen  so  that  each  beam  participates  equally  in  the 
vibration;  that  is,  if  you  assume  a  round  cross  section  of  beam,  the  maximum  stress 
level  experienced  in  beam  Li  under  gravity  loading  is  the  same  as  that  experienced 
in  L2. 

As  shown  in  Figures  4  and  5,  varying  the  length  of  the  first  and  second  links  does 
not  change  the  ratio  by  more  than  a  few  percent.  If  we  then  hold  the  lengths  constant 
and  change  the  weight  of  the  first  mass  and  the  stiffness  of  the  first  link,  as  shown  in 
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Lx 

L2 

Lg 

mi 

Kx 

1712 

A'2 

Link 

(cm) 

(cm) 

(cm) 

(N/cm) 

(kg) 

(N/cm) 

Steel 

45 

45 

HI 

12 

Aluminum 

60 

45 

■I 

|g 

80 

Table  1:  Parameter  Values  for  Beam  Comparisons 
Steel  Beams 


Actual 

to 

Estimated 

Frequency 


•X-  LI  =  23  cm 
LI  =  45  cm 
LI  =  68  cm 


Figure  4:  \’ariation  in  the  ratio  of  frequency  estimate  to  actual  natural  frequency  as 
a  function  of  link  length 

Figures  6  and  7,  the  change  is  much  larger  but  still  within  10%  of  the  original  ratio. 
Changing  the  first  ma.ss  has  the  largest  effect  on  the  estimate  but  even  then  fg  is  an 
estimate  good  to  within  20%  as  long  as  mi  is  kept  small. 

After  trying  different  values  of  lengths,  masses,  and  stiffnesses,  it  becomes  evident 
that  the  estimate  of  vibrational  frequency  based  on  static  gravitational  deflection  is 
very  good.  To  a  reasonable  approximation,  it  is  safe  to  say  that  the  vibrational 
frequency  estimate  based  on  endpoint  deflection  under  gravity  loading  is  good  to 
about  30%  of  the  actual  value.  As  fg  forms  a  lower  bound  to  the  vibrational  frequency, 
if  you  want  a  system  with  the  lowest  vibrational  frequency  for  a  given  endpoint 
deflection,  you  can’t  do  any  better  than  5/ y/Ff  (see  equation  7).  This  result  is  useful 
because  it  is  often  easier  to  calculate  an  endpoint  deflection  for  a  complex  system 
than  a  resonant  frequency. 

It  is  important  to  realize  that  the  gravitational  deflection  versus  natural  frequency 
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Aluminum  Beams 


Actual 

to 

Estimated 

Frequency 


•X'  LI  =  30  cm 
•o-  LI  =  60  cm 
L1  =  90  cm 


Figure  5:  Variation  in  the  ratio  of  frequency  estimate  to  actual  natural  frequency  as 
a  function  of  link  length 


Steel  Beams 


Actual 

to 

Estimated 

Frequency 


K1  (N/cm) 


X-  Ml  =  4  kg 
M1  =  7  kg 
♦-  M1  =  1 1  kg 


Figure  6:  Variation  in  the  ratio  of  frequency  estimate  to  actual  natural  frequency  as 
a  function  of  first  link  stiffness  and  joint  mass 
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Aluminum  Beams 


Actual 

to 

Estimated 

Frequency 


K1  (N/cm) 


•X-  M1  =  2  kg 
•o-  Ml  =  4.5  kg 
M1  =  7  kg 


Figure  7:  Variation  in  the  ratio  of  frequency  estimate  to  actual  natural  frequency  as 
a  function  of  first  link  stiffness  and  joint  mass 

relation  does  not  depend  on  the  length  of  the  beams  (see  equation  6).  If  you  want  a 
system  that  vibrates  at  1  Hertz,  you  must  live  with  a  endpoint  deflection  of  at  least 
25  centimeters  under  gravity,  regardless  of  how  long  your  beams  are.  Fortunately,  the 
required  endpoint  deflection  falls  off  quickly  as  the  vibrational  frequency  is  increased, 
(see  Figure  2)  The  robot  we  are  building  has  a  target  frequency  of  3Hz,  which  should 
give  it  a  little  over  a  one  inch  endpoint  deflection  due  to  gravity. 

Finally,  remember  that  this  relatir«pship  has  been  derived  based  on  the  endpoint 
deflection  under  a  gravitational  load.  Horizontal  estimates  are  found  just  as  easily 
if  you  assume  that  the  robot  has  been  placed  on  its  side  and  calculate  how  far  the 
endpoint  deflects. 

2.4  Length  of  the  Arms 

If  you  are  designing  a  robot  that  will  have  its  fundamental  frequency  of  vibration 
at  one  Hertz,  the  robot  is  going  to  sag  under  its  own  weight  at  by  least  25  cm.  A 
one  meter  arm  will  probably  not  work.  But  if  you  want  to  have  the  end  of  the  robot 
vibrating  with  an  amplitude  of  5  centimeters,  how  flexible  should  you  make  the  links? 
VVe  can  calculate  the  level  of  stress  inside  of  a  link  for  a  given  endpoint  deflection. 

In  terms  of  a  single  mass  system  with  a  constant  cross-section  beam  with  a  given 
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l  ij)  detlcciion  \vc  write  the  equatio 


a  = 


/ 


(Ifi) 


where  M  is  the  bending  moment  in  the  bar,  y  is  the  maximum  distance  from  the 
lu'utral  axis  and  /  is  the  moment  of  inertia.  The  bending  moment  is 


M  =  FL  =  K6u^L  =  —6u,  (IT) 

where  F  is  tiie  force  required  to  produce  an  end  deflection  and  L  is  the  length  of 
tl'.e  beam.  Combining  equations  (16)  and  (17),  we  have 


^  max 


^Ey 

L2 


(  IS) 


Equation  (18)  shows  that  the  stress  level  in  the  link  is  directly  proportional  to  the 
tip  deflection.  To  keep  the  bending  stresses  to  a  minimum  for  a  given  tip  deflection 
and  natural  frequency,  you  should  make  the  link  as  long  as  possible. 


3  Link  Types 

The  first  step  in  designing  a  flexible  robot  normally  includes  deciding  how  long 
the  links  will  be,  how  much  torque  will  be  available  at  the  joints,  and  what  the 
approximate  weights  and  sizes  of  the  joints  will  be.  Once  these  numbers  have  been 
selected,  you  need  to  pick  the  actual  shape  of  the  links  and  material  used  in  the  links. 
There  is  a  tradeoff  here  between  the  three  key  parameters:  stiffness,  mass,  and  stress 
level.  Typically  the  stiffness  will  have  been  set  by  the  decision  of  what  the  lowest 
frequency  of  vibration  will  be.  Then  you  must  design  links  for  this  stiffness  while 
keeping  the  stress  level  and  mass  of  the  link  to  a  minimum. 

3.1  Springs 

In  the  case  of  the  robot  we  are  currently  building,  it  worked  out  that  in  the  initial 
design  that  a  -3  Hertz  vibration  mode  put  uncomfortably  high  stresses  on  the  links. 
.At  the  time  we  were  calculating  based  on  links  made  out  of  solid,  straight  bars  of 
aluminum.  One  option  that  we  looked  at  to  eliminate  this  undesirable  stress  level  was 
to  replace  the  aluminum  bar  with  a  coil  spring.  At  first  glance,  the  coil  spring  seemed 
to  be  the  perfect  flexible  link;  not  very  stiff  and  able  to  undergo  large  deformations 
without  yielding. 

It  turns  out  that  springs  are  too  flexible  for  their  weight.  We  can  compare  the 
weight  and  flexibility  of  a  straight  link  made  out  of  a  bar  of  metal  with  diameter  d 
to  a  spring  made  out  of  the  same  bar  with  wire  diameter  d,  but  coiled  into  a  helix. 
Consider  the  overall  lengths  of  the  two  links  to  be  the  same. 
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First,  we  compare  the  bending  stiffness  of  the  spring  to  the  bending  stiffness  of 
the  bar.  The  general  equation  of  bending  for  a  spring  is 


r  Ed^p  I  1  \ 

0  ~  32LD  \  l  +  E/2g) 


(19) 


\\  here  r  is  the  bending  moment,  0  is  the  angle  of  deformation,  p  is  the  pitch  of  the 
spring,  L  is  the  length,  D  is  the  pitch  diameter  of  the  spring  and  E  and  G  are  the 
:i.  dulii  of  elasticity  and  rigidity.  If  we  assume  small  angles  of  deformation,  so  that 
=  dtj/dx  and  we  assume  that  the  spring  is  fi.xed  at  one  end  and  is  loaded  at  the 
c)rher  end  by  a  force  F  perpendicular  to  the  length  of  the  spring,  then  we  write  the 
deOection  of  the  end  of  the  spring  as  6,  where 

rL  T  E  f^- 

S  =  —  =  ~  (L  —  x)  dx  (20) 

Jo  K  K  Jo 


or 


(21) 


.Vote  ♦t'Ht  d  is  not  the  axial  deflection  of  the  spring  but  a  measure  of  how  far  the 
■'pi'ing  l)ends. 

.\ow  we  find  the  effective  spring  constant  of  the  spring  as 


K 


F 

S 


Ed*p 
32  PD 


(22) 


wlu-ie  we  have  made  the  approximation  E  —  2G.  For  the  round  bar  with  the  same 
diameter  d,  we  write 


3E/  37rEd^ 
'U'  ~  64E3 


(23) 


(  ’ombinining  equations  (22)  and  (23),  we  find  that 


_  '2p 

h'bar  SttD 


(24) 


I'or  a  practical  spring,  it  is  safe  to  assume  that  p  <  D,  so  K,  is  at  most  20%  of  Af.ar- 
The  weight  of  the  spring  can  be  compared  to  the  weight  of  the  straight  bar  .  The 
weight  of  the  spring  is  given  by 


W,  =  pxd^L 


and  the  weight  of  the  straight  bar  is 


W  bar  —  pTC d  Ij 


(25) 


(26) 
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(27) 


Then,  if  we  assume  that  (-n-Dfp)^  »  1,  we  get 

W,  _  ^ 

P 

Again,  it  is  safe  to  assume  that  p  <  Z),  so  the  weight  of  a  spring  made  of  a  coil  of 
wire  is  at  least  3  times  the  weight  of  a  link  made  of  a  straight  piece  of  that  w’ire. 

We  conclude  that  the  for  a  given  cross-section  d,  the  spring  weighs  at  least  three 
times  as  much  as  the  straight  bar  and  has  at  best  20%  of  the  bending  stiffness.  If 
you  want  a  very  low  bending  stiffness  and  don’t  care  about  weight,  the  spring  is  the 
way  to  go.  But  if  you  w^ant  a  given  stiffness,  then  the  spring  is  going  to  be  at  least 
15  times  as  heavy  as  a  straight  bar  of  metal.  Even  though  this  allows  you  to  not 
worry  about  breaking  your  link,  the  additional  weight  penalty  in  most  systems  is 
prohibitive. 


3.2  Material  Choice 


For  a  given  bending  stiffness,  the  material  used  for  the  link  will  determine  its 
size,  weight,  and  how  much  deflection  it  will  undergo  before  yielding.  A  material 
with  a  high  yield  stress  is  not  necessarily  the  best  choice.  If  it  has  a  high  modulus 
of  elasticity,  the  link  will  have  to  be  thinner  to  get  the  same  bending  stiffness.  The 
smaller  moment  of  inertia  will  result  in  a  higher  stress  level  than  that  of  a  link  made 
of  a  material  with  a  high  yield  strength  and  a  low  modulus  of  elasticity.  A  similar 
situation  exists  with  weight. 

To  compare  materials,  we  first  a.ssume  that  the  link  is  a  solid  bar  of  metal  with  a 
round  cross  section  of  diameter  d.  We  choose  a  round  cross  section  for  two  reasons: 
First,  it  has  the  same  stiffness  when  bent  in  any  direction,  where  the  stiffness  is 
deflection  of  the  end  of  the  link  with  respect  to  the  force  applied  at  the  end.  Second, 
for  a  given  bending  moment  and  stiffness  K,  a  solid  bar  experiences  a  lower  bending 
stress  than  a  hollow  bar.  (Consider  two  bars  with  the  same  moment  of  inertia;  one 
solid  and  one  hollow.  The  hollow  bar  will  have  a  larger  diameter  and  the  bending 
stress  is  directly  proportional  to  the  diameter.)  If  we  assume  that  we  know  the 
length  of  the  link,  the  stiffness,  and  the  applied  bending  moment,  we  can  calculate 
the  stress  and  weight  of  the  link  as  a  function  of  density,  modulus  of  elasticity  and 
yield  strength. 

Start  with  the  standard  formulas 


(28) 


and 


(29) 
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where  E  is  the  modulus  of  elasticity  of  the  metal,  M  is  the  applied  bending  moment, 
/  is  the  moment  of  inertia  of  a  round  cross  section  and  L  is  the  length  of  the  beam. 
Rewrite  equation  (29)  as 


(3U) 


or 


d  = 


64  £3  A' 
37rA 


(31) 


Fo  find  the  stress  in  the  beam,  combine  equations  (28)  and  (31) 


6M  .IKE^ 
KL^\I  UttL 


(32) 


Ec^uation  (32)  shows  that  with  £,  K  and  M  fixed  for  a  beam,  the  material  with  the 
lowest  fraction  of  stress  to  yield  stress  will  be  the  one  with  the  largest  value  of  V’, 
where  xl'  is  defined  as: 


^  = 


^ yield 


(33) 


The  larger  the  value  of  the  larger  the  bending  moment  the  link  will  withstand 
before  yielding. 

Another  material  comparison  is  the  weight  of  the  bar  as  a  function  of  the  length 
and  the  spring  constant.  Using  equation  (31)  and  taking  p  as  the  density  of  the  link, 
we  get 


npxPL  vL^K 

=  — 


(34) 


where  W  is  the  weight  of  the  link.  For  a  given  L  and  A',  the  beam  with  the  lowest 
wf'ight  will  be  the  one  that  minimizes  <f>,  where  is  given  by: 


(35) 


1  he  larger  the  value  of  <i>,  the  heavier  link. 

The  optimal  link  material  would  have  a  large  xjj  and  a  small  d>.  In  practice,  no 
material  is  optimal  although  Titanium  comes  close.  Table  2  shows  the  value  of  V’ 
and  (f)  worked  out  for  a  number  of  different  common  metals.  Titanium  has  the  best 
strength  to  stiffness  and  2014-T6  Aluminum  is  a  practical,  lightweight  alternative 
although  its  resistance  to  fatigue  is  limited. 


4  Joint  Flexibility 

Robots  have  two  types  of  flexibility;  in  the  joints  and  in  the  links.  Joint  flexibility 
ap])ears  to  the  system  as  springs  in  series  with  the  links.  To  illustrate  this,  look  at 
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Material 

(MPa) 

E 

(GPa) 

(kg/m^) 

1100-0  Aluminum 

34 

1.4 

330 

2014-T6  Aluminum 

415 

16.8 

330 

6061-T6  Aluminum 

275 

72 

2800 

11.1 

330 

1015  Steel 

324 

207 

■■■ 

5.9 

535 

4140  Steel 

655 

207 

12.0 

535 

Magnesium  Alloy 

240 

45 

■■ 

13.8 

268 

Titanium  Alloy 

830 

114 

4400 

23.4 

412 

Table  2:  Link  Material  Comparison 

the  two  link  model  shown  in  Figure  3.  Now  add  in  a  torsional  spring  of  value  Ki 
where  the  first  link  meets  the  wall  and  a  torsional  spring  of  value  K2  between  the 
links.  Equations  (9),  (10),  and  (11)  relate  the  deflection  of  a  point  to  a  unit  force 
applied  at  a  point,  so  we  can  modify  them  to  take  into  account  the  torsional  springs: 


/  1 

All  —  On  + 

(36) 

_ / _ .  +  -^i/2 

a  12  —  a  21  0,21  + 

(37) 

,  Li  +  L2  +  Lj  _  Z-2  +  Tj/2 

—  <^22  1  + 

K2 

(38) 

These  modified  flexibility  values  can  be  substituted  in  equations  (13)  and  (15)  to 
give  the  estimates  for  the  lowest  natural  frequency  of  the  system.  We  still  have  fg 
as  a  lower  bound  for  /„,  and  assuming  that  the  springs  are  not  extremely  flexible,  fg 
still  gives  a  good  approximation  of  the  lowest  natural  frequency. 

5  Conclusion 

The  endpoint  deflection  of  a  robot  in  a  gravity  field  due  to  bending  in  the  links 
and  flexiblity  in  the  joints  forms  a  lower  bound  for  the  natural  frequency  of  vibration 
of  the  system.  When  all  of  the  links  are  flexible  and  the  weight  of  the  payload  is  not 
negligible  in  comparison  to  the  weight  of  the  joints,  this  estimate  is  very  close  to  the 
true  frequency.  As  this  estimate  forms  a  lower  bound,  if  your  goal  is  to  minimize 
the  fundamental  natural  frequency  of  the  system,  you  should  make  links  that  are 
as  long  aa  possible  and  have  small  masses  at  the  joints.  The  “tip  deflection  under 
gravity”  calculation  is  not  restricted  to  the  vertical  case;  we  can  put  an  imaginary 
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liiavitational  field  pointing  in  any  direction,  calculate  the  deflection  of  the  robot,  and 
<'stiniate  the  natural  frequency  of  vibration. 

For  maximum  flexibility  with  minimum  stiffness,  the  best  material  choice  for  the 
links  is  titanium.  Some  aluminum  alloys  may  be  inexpensive  alternatives,  but  the 
fatigue  characteristics  of  aluminum  are  unacceptable  in  most  applications.  Coiled 
springs  give  very  little  stiffness  in  comparison  to  their  strength,  but  their  weight 
makes  them  unusable  in  all  but  planar  applications. 

There  are  physical  difficulties  to  building  a  test  fixture  with  a  very  low  frequency 
of  vibration.  To  keep  stress  levels  low  in  a  flexible  robot,  long  links  are  prefered.  but 
may  not  be  possible  if  the  robot  must  fit  inside  of  a  laboratory.  The  longer  links 
have  a  larger  moment  of  inertia  and  need  larger  motors  to  drive  them  around.  Joints 
that  contain  motors  carry  a  great  deal  of  weight  in  the  form  of  motors,  gears,  brakes, 
encoders,  and  ball  bearings.  Trying  to  optimize  the  robot  to  have  low  joint  masses, 
long  links,  and  adequate  torque  at  the  joints  is  difficult.  With  these  constraints  in 
mind,  we  are  currently  assembling  a  three  degree  of  freedom  flexible  robot.  It  has 
a  lowest  natural  frequency  of  vibration  of  under  3  Hertz,  a  payload  capacity  of  3 
pounds  and  a  total  reach  of  at  least  51  inches.  We  expect  that  it  will  be  operational 
by  June  of  1989. 

References 

[1]  Cannon,  Robert  H.  and  Schmitz,  Eric, 

“Initial  Experiments  on  the  End-Point  Control  of  a  Flexible  One-Link  Robot”, 
The  International  Journal  of  Robotics  Research,  Vol.  3,  No.  3.,  Fall  1984. 

[2]  Alberts,  T.E.,  Hastings,  G.G.,  Book,  W.J.,  and  Dickerson,  S.L., 

“Experiments  in  Optimal  Control  of  a  Flexible  Arm  with  Passive  Damping”, 
Fifth  VPISSU/AIAA  Symposium  on  the  Dynamics  and  Control  of  Large  Struc¬ 
tures,  1985. 

[3]  Hollars,  Michael  G.  and  Cannon,  Robert  H,, 

"Initial  Experiments  on  the  End-Point  Control  of  a  Two-Link  Manipulator  with 
Flexible  Tendons”,  ASME  Winter  Annual  Meeting,  November  19,  1985. 

[4]  Tilley,  S.W.,  Cannon,  R.H.,  and  Kraft,  Ray, 

“End  Point  Force  Control  of  a  Very  Flexible  Manipulator  with  a  Fast  End 
Effector”,  ASME  Winter  Annual  Meeting,  December,  1986. 

[5]  Bayo,  E., 

“Computed  Torque  for  the  Position  Control  of  Open-Chain  Flexible  Robots”, 
IEEE  International  Conference  on  Robotics  and  Automation,  April  25  29, 
1988. 


14 


[6]  Pfeiffer,  F.  and  Gebler,  B., 

“A  Multistage- Approach  to  the  Dynamics  and  Control  of  Elcistic  Robots”, 
IEEE  International  Conference  on  Robotics  and  Automation,  April  25-29, 
1988. 

[7]  Blevins,  Robert  D., 

Formulas  for  Natural  Frequency  and  Mode  Shape,  Van  Nostrand  Reinhold 
Company,  Appendix  A,  1979. 

[8]  Thomson,  William  T., 

Theory  of  Vibration  with  Applications,  Prentice-Hall,  Inc.,  1981. 


15 


DISTRIBUTION: 


Defense  Technical  Information  Center 

Computer  Sciences  Division 
ONR,  Code  1133 

Navy  Center  for  Applied  Research  in  Artificial  Intelligence 
Naval  Research  Laboratory,  Code  5510 

Dr.  A.L.  Slafkosky 

Headquarters,  U.S.  Marine  Corps  (RD-1) 

Psychological  Sciences  Division 
ONR,  Code  1142PT 

Applied  Research  &  Technology 
ONR,  Code  12 

Dept,  of  the  Navy 
Naval  Sea  Systems  Command 
NAVSEA  90 

Dr.  Charles  Schoman 

David  Taylor  Naval  Ship  R&D  Center 

NSRDC  18* 


